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Abstract 

The temperature and polarization inhomogeneities of the Cosmic Microwave Background 
might bear the mark of pre-decouphng magnetism. The parameters of a putative magnetized 
background are hereby estimated from the observed temperature autocorrelation as well as 
from the measured temperature-polarization cross-correlation. 
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In recent years the observations of the Cosmic Microwave Background (CMB in what 
follows) have been translated into a set of constraints on the standard parameters of an 
underlying ACDM scenario where A stands for the dark energy component (parametrized 
in terms of a cosmological constant) and CDM is the acronym of cold dark matter which 
is the dominant source of non-relativistic matter of the model. More specifically, the five 
year release of the WMAP data [T] (WMAP 5yr in what follows) can be used to infer the 
parameters of a putative ACDM paradigm where flxo = Pxo/Pcrit denotes the (present) 
energy density of the species X in units of the critical energy density pcrit- The best-fit 
values of the Qxo as well as of the other parameters of the model can be determined, for 
instance, from the WMAP Syr data alone □ (see also [2] for earlier results) and they are 

f^bo = 0.0441 ± 0.0030, fico = 0.214 ±0.027, ^ao = 0.742 ± 0.030, (1) 

where the subscripts b and c stand, respectively, for baryons and CDM components; Eq. ([T]) 
does not exhaust the list of parameters which also include the Hubble constant (i.e. -f^o)) 
the optical depth of reionization (i.e. ere) and the spectral index of (adiabatic) curvature 
perturbations (i.e. rig): 

Ho = 71.9^^ 7 km/sec Mpc'S e^e = 0.087 ± 0.017, = 0.963'^°^o\5; (2) 

ere corresponds, in the framework of the WMAP Syr data, to a typical redshift Zj-e = 11 ±1.4. 
In the ACDM paradigm the temperature and polarization inhomogeneities are solely sourced 
by the fluctuations of the spatial curvature 71 whose Fourier modes obey 

-» 27r^ -> / k \ "s-i 

{n{k)n{p)) = ^6^'\k + p), Vn{k) = An[^J , (3) 

where An = (2.41 ± 0.11) x 10~^ is the amplitude of the curvature power spectrum at the 
pivot scale kp = 0.002 Mpc"^ The success of the ACDM paradigm suggests the possibility 
of constraining also other (non-standard) parameters such as a background of relic gravitons 
[Ij, non-adiabatic modes possibly present prior to equality [5] or even a background of relic 
magnetic fields whose presumed existence is not in contrast with the actual observations 
of large-scale magnetism in galaxies |6], clusters [7] or even (with mandatory caveats) in 
superclusters [8]. Can we falsify, at a given confidence level, the hypothesis that large-scale 
magnetic fields are already present prior to decoupling? 

This paper follows a series of investigations ^ where large-scale magnetic fields have 
been consistently included at all the stages of the Boltzmann hierarchy and within a faith- 
ful plasma description incorporating the evolution of the space-time curvature and of its 
inhomogeneities. For frequencies smaller than the (electron) plasma frequency, the Ohmic 
current is solenoidal and the electric fields are suppressed by (Coulomb) conductivity: in 

^The combination of the WMAP 5yr data with other cosmological data sets (e.g. supernovae |3] and 
large scale structure 4 observations) lead to compatible values pj. 
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this regime the effective dynamical variable is the centre of mass velocity of the electron ion 
fluid (sometimes named, with confusing terminology, the baryon velocity). For frequencies 
larger than the plasma frequency, the propagation of electromagnetic disturbances follows a 
two-fluid approach where ions and electrons are treated as independent dynamical entities: 
Faraday effect induces, in this regime, a B-mode polarization [10] which might be relevant 
when (and if) compelling experimental determinations of the B-mode autocorrelation will 
be available. 

The temperature autocorrelations (TT spectra in what follows) and the temperature- 
polarization cross-correlations (TE spectra in what follows) are determined with reasonable 
accuracy but, so far, there has not been any attempt to extract the magnetic field parameters 
from the TT and TE angular power spectra. The modest purpose of the present investigation 
is to fill such a gap by using the WMAP 5yr data release in the light of the minimal version of 
the magnetized ACDM scenario (mACDM in what follows) where the initial conditions of the 
Einstein-Boltzmann hierarchy are assigned in terms of a magnetized adiabatic mode [9] and 
the (magnetized) TT and TE correlations are determined accordingly. The characteristic 
parameters of the (minimal) mACDM scenario are the magnetic spectral index tt-b and the 
(comoving) magnetic field intensity i?L whose associated energy density can be referred to 
the photon background as Qbl = B^/ i^'^P-y) ■ Dividing ^bl by An we get 



where T^o = 2.725±0.001 K is the temperature of the CMB and the comoving magnetic field 
i?L is measured in nG (i.e. InG = 10^^ Gauss = 6.924 x lO^^^GeV^). Since the magnetic 
fields are stochastically distributed, the ensemble average of their Fourier modes obeys: 

{B,{k)Bj{p)) = —P,^{k)PB{k)5^'\k + p), VB{k) = Ab[-j^) , (5) 

where Pij{k) = {k'^6ij — kikj)/k'^; Ab the spectral amplitude of the magnetic field at the 
pivot scale = Mpc~^ [9l [10]. In the case when ub > 1 (i.e. blue magnetic field spectra), 
Ab = (27r)"B-i5^/r[(nB - l)/2]; if tzb < 1 (i.e. red magnetic field spectra), ^b = [(1 - 
Ub) /2]{kA/ ki^)^^'"^^^ Bl where /ca is the infra-red cut-off of the spectrum. In the case of white 
spectra (i.e. ub = 1) the two-point function is logarithmically divergent in real space and 
this is fully analog to what happens in Eq. Q when Ug = 1, i.e. the Harrison- Zeldovich 
(scale-invariant) spectrum. By selecting k^^^ of the order of the Mpc scale the comoving 
field i?L represents the (frozen-in) magnetic field intensity at the onset of the gravitational 
collapse of the protogalaxy [S]. 

The evolution and regularization of the magnetic fields involve an ultraviolet cut-off k-o 
(related to the thermal diffusivity scale at last-scattering) as well as an infra-red cut-off /ca 
related to the (comoving) angular diameter distance at last scattering. The values of kA and 
/cd do depend upon the best-fit parameters of Eqs. ([I])-([3|). The redshift of last scattering 
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2*, for typical values of the parameters close to the values of Eqs. ([I])-®, turns out to be: 

. -fonr 00. 0.0783 c^r°-2^^ 0.560 , , 

= 104811 + A., ^.= |i^3a,5(,>„l . 1 + 21.1.... . (6) 

where um = uji, + Uc and u\y = /io^bo, = ^o^co; in Eq- dS]) /i = 1.24 x 10"^ and 
/2 = 0.738. In the case of the parameters of Eqs. ([l])-(l3]), Eq. ([6]) implies 2;* = 1090.5 in ex- 
cellent agreement with the estimate of the WMAP collaboration P], i.e. z^, = 1090.51 ± 
0.95. From the value of it is possible to compute tr^, = Pr(2;*)/pm(-2*) = 4.15 x 
10^^/c<jm (-2*/10'^). Recalling that the comoving angular diameter distance can be written 
as Da{z^) = 2dA{z^) / (Hq^/Qmo) , kx{z^) and /i;d(-2*) are determined in terms of Eqs. ([T])-(j3]) 

k^{z,) 2240rfA(2;*) f z,\^l^ -0.11 



where, for the best fit parameters hsted in Eqs. ([I])-®, dx{z^) = 0.8569 leading to Da{z^) = 
14110.8 Mpc in excellent agreement with the WMAP 5yr determination (i.e. Da(-2*) = 
141151.^1891 Mpc) implying kA{z^) = 7.08 x 10^^ Mpc"\ 

Throughout the numerical scan of the parameter space the initial conditions of the Boltz- 
mann hierarchy have been given in terms of the magnetized adiabatic mode which has been 
discussed in [9l [10] . The numerical code employed here is an improved version of the one used 
in [9l [To] and it is based on Cmbfast [11] (which follows, in turn, Cosmics [12]). Let us now 
remark that the joined two-dimensional marginalized contours for the various cosmological 
parameters identified already by the analyses of the WMAP 3yr data [2] are ellispses with an 
approximate Gaussian dependence on the confidence level as they must be in the Gaussian 
approximation (see, e.g., already the first paper of Ref. [2] and, in particular. Fig. 10). In 
the light of this observation it is then plausible to determine the confidence intervals for the 
2 supplementary parameters of the model (i.e. n-g, Bi^) by using an appropriate gridding 
approach while the remaining parameters of Eqs. are assumed to be known and 

follow a Gaussian probability density function. This approach has been also followed when 
discussing, some time ago, the impact of non-adiabatic modes on the WMAP 3yr data (see, 
e. g., second paper of Ref. [5]). In the numerical fit to the data the shape of the likelihood 
function can be determined by evaluating the least square estimator 



X^(nB,i?L) = E 



d°^^)_Q(nB,i?L) 



(obs) 
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(8) 



i 

where a"^^ are the estimated errors from the observations for each multipole and where the 
functional dependence of Ci{nB, Bi) is given by the underlying theory (i.e. the magnetized 
ACDM model) which we try to falsify by comparing its predictions to the actual observations. 
The observed angular power spectra (i.e. C^^^) are derived by using the (further) estimators 
and i.e. 



C, 



(TT) _ 1 ^ |.(T)|2 ^(TE) _ 1 ^p |.(T).(E)*| 
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whose distribution becomes Gaussian, according to the central hmit theorem, when £ — > oo. 
The true a~^^ and a'-^-' are not accessible to direct observation, however, their estimates 
(i.e. a^"J2 and ci^^) can be directly inferred from the measured temperature and polarization 
inhomogeneities [1]. The minimization of Eq. ([8]) is equivalent to the minimization of the 
lognormal likelihood function C = —2 InL where L is given by 

L(data| ne, B^) = L^^^e->^"^^^^^^^^/\ (10) 

Thus, the minimization of Eq. ([8]) is equivalent to the maximization of the likelihood of Eq. 
f fTOj) . In Fig. [T]the filled (ellipsoidal) spots in both plots highlight the minimal value of the 



nB= 1.598, BL = 3.156nG nB= 1.616, Bl = 3.218 nG 
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Figure 1: Likelihood contours in the plane (ne, -Bl)- In both plots the inner curve cor- 
responds to Ax^ = 2.3 (i.e. 68.3% or la of likelihood content) while the outer curve 
corresponds to A^^ = 6.17 (i.e. 95.4% or 2 a of likelihood content). 

(i-e. Xmin) while the contour plots represent the likelihood contours in the two parameters 
n-Q and -Bl- In both plots of Fig. [T]the boundaries of the two regions contain 68.3% and 
95.4 % of likelihood as the values for which the has increased, respectively, by an amount 
Ax^ = 2.3 and Ax^ = 6.17. In the plot at the left of Fig. [T]the data points correspond to 
the measured TT correlations contemplating A^tt = 999 experimental points from ^ = 2 to 
£ = 1000. In the plot at the right the data points include, both, the TT and TE correlations 
and the total number of data points increases to N^y + N'^-e = 1998. The values of ne and -Bl 
minimizing the when only the TT correlations are considered turn out to be tib = 1.598 
and -Bl = 3.156 nG (see Fig. [T]plot at the left); in this case the reduced is 1.09. When also 
the TE correlations are included in the analysis the reduced x^ diminishes from 1.09 to 1.03 
and the values of ub and -Bl minimizing the x^ become ub = 1.616 and -Bl = 3.218 nG. In 
the frequentist approach, the boundaries of the confidence regions represent exclusion plots 
at 68.3 % and 95.4 % confidence level. In this sense the regions beyond the outer curves in 
Fig. [Hare excluded to 95% confidence level. 
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In Fig. [2] the full line corresponds to the WMAP Syr best fit of Eqs. (P)-©. The 
(hardly distinguishable) dashed lines, in both plots, illustrate the niACDM result for {n-Q = 
1.616, Sl = 3.218 nG) corresponding to Xmin; determined from the joined analysis of the 
TT and TE angular power spectra (see plot at the right in Fig. [1]). In both plots of Fig. [1] 
the (binned) data points are reported §. The impact of the magnetic field parameters can 
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Figure 2: The TT and TE angular power spectra for the best-fit parameters stemming from 
the right plot of Fig. [T] (dashed line) and from the best fit to the WMAP Syr data alone 
(full line) . The data points of the WMAP Syr release are included in their binned form. 

be appreciated by monitoring the relative ratios of the first three peaks. Using the notation 
Qe = i{i + l)Ci/ (27r), the relative heights of the acoustic peaks computed in the case of the 
best-fit of Fig. [2] are: 

Hi = -^ = 6M2, H^ = ^ = 0AA7, H^ = ^ = 0.981, (11) 

where = 220, ^2 = 53S and £3 = 816 are, respectively, the locations of the first three acous- 
tic peaks. In the the case of Eqs. ([I])-([3]) (i.e. in the absence of a magnetized background) 
the same ratios are, respectively. Hi = 6.876, H2 = 0.447 and = 0.963. In Fig. [3] (plot 
at the left) the TT power spectra are illustrated with the same value of the spectral index as 
in Fig. [2] (i.e. ub = 1.616) but with a magnetic field intensity i?L = 6 nG which is more than 
two a away from Xmin (see Fig. [H plot at the right): as expected, rough inspection of the 
left plot in Fig. [3l reveals that the mACDM estimate does not correctly reproduce the data 
when (ne = 1.616, i?L = 6 nG) especially around the first and third peaks. In the latter case 
the the relative ratios of the peaks become {Hi, H2, H3) = (7.142, 0.444, 1.012): i.e. Hi and 

The unbinned data (which are the ones used in the analysis) contemplate all the multipoles from £ ~ 2 
to £ = 1000 both for the TT and for the TE (observed) power spectra. The binned data contain instead 34 
(effective) multipoles in the TE correlation and 43 (effective)multipoles in the TT spectrum. Following the 
usual habit, to make the plots more readable, the binned data points have been included in Fig. [2l 
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i^s exceed the best fit while H2 is comparatively lower than in Eq. flTT]) . This is a general 
trend in the excluded region. Consequently, the whole approach of ^9J is rather effective in 
constraining variations of i?L in the nG range. Similarly, by moving the magnetic spectral 
index in the excluded region of the parameter space the distortions of the peaks jeopardize 
the goodness of the fit. In Fig. [3] (plot at the right) the magnetic pivot wave-number has 
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Figure 3: Left plot: the TT power spectra which are more than 2a away from the best-fit 
point (full line) are compared with the best fit of Fig. [2] (dashed line). Right plot: contour 
plots in the (ne, -Bl) plane for /cl = 0.1 Mpc~^ (see also Fig. [T]where /cl = 1 Mpc~^). 

been moved from /cl = 1 Mpc~^ to /cl = 0.1 Mpc~^, consistently with the indetermination of 
the scale of protogalactic collapse. A reduction of the pivot wave-number implies that the 
magnetic field is regularized over a broader window in real space; hence the allowed ampli- 
tude of the putative magnetic field intensity is reduced. In the right plot of Fig. [3] the is 
minimized when n-Q = 1.617 (i.e. comparable with the value of ne deduced in Fig. [2] when 
ki^ = lMpc~^) and, as anticipated, for a smaller magnetic field, i.e. i?L = 1.578 nG. The 
modification in the pivot wave-number also modifies the shape of the confidence contours 
(compare, e.g.. Fig. [Hand Fig. [3]). 

The large-scale magnetic field possibly present prior to decoupling have been estimated, 
for the first time, from the temperature autocorrelations as well as from the cross-correlations 
between temperature and polarization in the case of the magnetized adiabatic mode [9]. The 
obtained results prove in explicit terms a guessed degeneracy between the spectral index 
and the magnetic field intensity. In a frequentistic perspective, beyond the outer contours 
of Figs. [1] and [3] the parameter space of the (minimal) mACDM scenario is excluded to 95% 
confidence level. Global observables (such as, for instance, the heights and depths of the 
acoustic oscillations) could be probably used as a normal parameter basis to infer an explicit 
analytic dependence of the TT and TE correlations upon the parameters of the magnetized 
background. In this direction work is in progress. 
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